Molecular Probe of Pairing in the BEC-BCS Crossover 
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We have used optical molecular spectroscopy to probe the many-body state of paired 6 Li atoms 
near a broad Feshbach resonance. The optical probe projects pairs of atoms onto a vibrational level 
of an excited molecule. The rate of excitation enables a precise measurement of the closed-channel 
contribution to the paired state. This contribution is found to be quite small, supporting the concept 
of universality for the description of broad Feshbach resonances. The dynamics of the excitation 
provide clear evidence for pairing across the BEC-BCS crossover, and into the weakly interacting 
BCS regime. 

PACS numbers: 03.75.Ss, 33.80.Ps, 67.90.+Z, 03.75.Nt 



Feshbach resonances were recently used to explore the 
possibility of superfluidity and Cooper pairing in atomic 
Fermi gases of 40 K and 6 Li. Two experiments have pro- 
duced evidence for correlated pairs while others re- 
port evidence for superfluidity Q. These experiments 
are noteworthy because they are performed on the high- 
field side of the resonance, where the s-wave scattering 
length a is negative and, according to two-body physics, 
diatomic molecules associated with the resonance are en- 
ergetically unstable. These intriguing experiments raise 
important questions regarding the nature of Feshbach- 
induced pairs, including their relation to bound molecules 
and the role played by many-body effects. 

Atomic Feshbach resonances involve a collision of a 
pair of atoms in an open channel coupled to a bound 
molecular state of a closed channel |3| . In the broad 6 Li 
resonance, the open channel corresponds to two atoms in- 
teracting mainly via the electronic spin-triplet potential, 
while the closed-channel is predominantly spin-singlet. 
Two-body theory predicts the formation of a weakly- 
bound molecular state on the "BEC" side resonance, 
where a > 0. These "dressed" molecules are superpo- 
sitions of open-channel atoms with the "bare" molecules 
of the closed channel. For resonances that are broad com- 
pared to the Fermi energy, the closed channel character of 
the dressed molecules is expected to be small throughout 
an experimentally relevant region about the resonance 
0. In this case, the resonance may be well-described 
by a single-channel model where the physics is universal, 
such that the macroscopic properties of the superfluid 
are independent of the microscopic physics that underlie 
the two-body interactions. Universality in the descrip- 
tion of pairing in atomic gases establishes their relevance 
to other systems, most notably high-temperature super- 
conductors. 

In the experiment reported here, a laser is used to 
project the dressed molecules/pairs onto an excited sin- 
glet molecular state. By starting with an evaporatively 
cooled gas on the BEC side of the resonance, followed by 
an adiabatic change in the magnetic field, a nearly zero 
temperature gas can be probed throughout the BEC- 



BCS crossover. This enables a direct measurement of 
the closed channel fraction and, for the first time in an 
atomic gas, provides clear evidence for the presence of 
pair correlations in the weakly interacting BCS regime. 

Our apparatus and methods for producing strongly in- 
teracting Fermi gases of 6 Li atoms have been described 
previously [f|. The atoms are confined in a single-beam 
optical trap which, at full laser intensity, has a trap depth 
of 25 fiK and radial and axial freqencies of v r = 2270 Hz 
and v z = 21 Hz, respectively. Curvature in the mag- 
netic bias field modifies the axial frequency such that 
v' z = y/uf+XB, where A = .029(3) Hz 2 /G. An inco- 
herent equal mixture of the two lowest going Zeeman 
sublevels (F = 1/2, tuf = ±1/2) is created |5J. For the 
magnetic fields of interest (B > 600 G) these states are 
nearly electronically spin polarized, but differ in their nu- 
clear spin projections. They interact strongly via a broad 
Feshbach resonance located near 834 G |6J, [7| • 

A molecular BEC is created as in other 6 Li experiments 
. We prepare the spin mixture at 754 G (corresponding 
to a ~ 3680 a a and dressed molecules form through 
three-body recombination. The molecules are evapora- 
tively cooled by reducing the trap depth in 750 ms in an 
approximately exponential trajectory. Their distribution 
is imaged in the trap with the same optical probe used 
to image the atoms. Figure shows that the condensate 
fraction can be controlled by the final trap depth. The 
crossover regime is probed by preparing the gas at 754 G 
and adiabatically changing the magnetic field to any de- 
sired final value. We find that adiabatic field excursions, 
even those going across the resonance, result in no de- 
tectable heating upon returning to the original field [ljj • 

A molecular probe laser drives transitions between 
the dressed molecules/pairs and an electronically excited 
molecular singlet level. These transitions result in spon- 
taneous emission and a detectable loss of atoms from the 
trap. We chose to drive transitions to the v' = 68 level of 
the A 1 E+ excited state because it has the largest Franck- 
Condon wavefunction overlap (0.077) with the resonant 
channel bare molecule, the X 1 S+(z; = 38) level. The 
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FIG. 1: In-situ absorption image profiles showing a molecular 
BEC. These images were recorded at a field of 695 G after 
evaporation at 754 G. For (a), the optical trap depth was 
lowered to 0.5 fiK. The solid line is a fit to a Gaussian (dotted 
line) plus Thomas-Fermi distribution which distinguishes the 
condensate from residual thermal molecules. For (b), the trap 
depth is reduced to 0.27 /J.K, producing an essentially pure 
molecular condensate with the number of molecules, N = 
46, 000. We estimate that the condensate fraction is >90%, 
implying T/T c < 0.5, where T c is the critical temperature for 
BEC at 695 G. The solid line is a Thomas-Fermi distribution. 



v' — 68 level has a classical turning point of 36 clq. A 
quantity Z is denned in terms of T, the rate of pho- 
toexcitation, by T — Zil 2 /j, where fl = (i/j v '—es(S = 
0)\d ■ E L \tp v=3S (S = 0)) is the on-resonance Rabi fre- 
quency, d is the transition dipole, E^ is the laser field 
of the molecular probe, and 7 = (2tt) 11.7 MHz is the 
linewidth of the excited molecular state 01 • The dressed 
molecules/pairs can be expressed as a superposition of 
the v = 38 singlet molecules and free atom pairs in the 
triplet channel [3j: 

IW> = ^ 1/2 |^= 38 (5 = Q)) + (l-Zyl 2 \US = 1)), (1) 

and Z can be identified as the component fraction of the 
dressed molecules in the closed-channel. For a typical 
probe laser intensity of 20 mW/cm 2 , Q = (2tt) 2.6 MHz 

Z is measured at various fields by the following process: 
first, the molecular gas is prepared at 754 G. The field 
is then linearly ramped to a new value B in 160 ms and 
held for 20 ms before the molecular probe is pulsed on 
for a fixed duration. The probe removes atoms at a rate 
r determined by Z . The field is then linearly returned to 
754 G in 160 ms, where the number of remaining atoms is 
determined by optical absorption imaging. The measured 
numbers are normalized to data obtained by the same 
procedure, but without firing the molecular probe. 

Figure|21shows the normalized loss of atoms at B = 695 



FIG. 2: Loss of signal vs. molecular probe duration at 695 
G. The open circles correspond to a gas evaporatively cooled 
to a nearly pure molecular BEC, while the closed circles cor- 
respond to full trap depth, where T/Tf — 0.75. The dashed 
lines are fits to exponentials, with a leftover fraction of 25% 
in the high temperature case. The time axis for the BEC data 
was scaled to account for differences in molecular probe laser 
intensity between the BEC data (30 mW/cm 2 ) and the high 
temperature data (15 mW/cm 2 ). We have verified that the 
loss rate depends linearly on intensity. The error bars repre- 
sent the statistical standard deviation from the mean of ~10 
independent measurements. 



G as a function of the molecular probe duration for 
two different temperatures. In the case of the pure 
condensate, the gas is well into the BEC regime, since 
a ~ 1510 a Q and kpa ~ 0.15. As expected in this case, 
the entire trap can be depleted for sufficiently long probe 
duration. Furthermore, the signal decays exponentially, 
indicative of a one-body loss process. For the higher tem- 
perature data, the loss is initially exponential, but 25% 
of the initial number remain after a long probe duration. 
The remainder can be understood as the dissociated atom 
fraction, which is determined by the fitted temperature 
T = 1.8 fiK and the calculated binding energy E\, = 13 
fiK of the dressed molecules. A detailed balance calcula- 
tion gives approximate agreement with the observed dis- 
sociated fraction. The probe only weakly couples to free 
atoms since they have a relatively small singlet character 
and excitation occurs only by two-body photoassociation. 

Figure shows the loss of signal vs. probe duration 
at 865 G, where a ~ —15600 a . At this field, the gas is 
in the strongly interacting regime where kp\a\ > 1. Ac- 
cording to two-body physics, there are no bound states 
above the resonance at 834 G, and probe-induced loss 
would arise exclusively from two-body photoassociation. 
Nonetheless, the observed loss for the gas prepared in a 
nearly pure BEC fits to an exponential rather than to a 
two-body process. At full trap depth, the decay consists 
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FIG. 3: Same as for Fig. H except at 865 G. The dashed 
line in the case of the full trap depth data (closed circles) is 
a fit to a "two-fluid" model where one component decays via 
a rapid one-body loss process and the other via a slower two- 
body loss process. Approximately 75% of the gas is lost by 
the initial fast process. 



of two parts: an initial exponential decay, followed by a 
much slower two-body process. As the temperature of 
the cloud is T ~ 0.75 Tp, this initial fast decay may in- 
dicate the presence of uncondensed paired fermions 0] 
or finite-lifetime molecules. The slower (two-body) pro- 
cess is ascribed to free-bound photoassociation, which is 
supported by the fact that the extracted two-body rate 
coefficient, K 2 = 4.9(3.3) x 1CT 10 (cm 3 s" 1 )/(W cm" 2 ), 
agrees well with the calculated value of 9.8(2.6) x 10~ 10 
(cm 3 s _1 )/(W cm -2 ) obtained using the expression for 
K2 given in Ref. , where the uncertainties arise mainly 
from the temperature determination. 

Figure 0| shows the extracted values of Z for fields be- 
tween 600 and 920 G. Below 600 G, the dressed molecule 
lifetime is too short to obtain a reliable measurement of 
Z . Also shown in Figure 01 are the results of a coupled 
channels calculation, which represents an exact two-body 
theory. An analytic expression for Z on the BEC side of 
the resonance has been given in Ref. and is in good 
agreement with our calculation. While the two-body the- 
ory accurately represents the data for fields below reso- 
nance, there is complete disagreement above resonance. 
Two-body theory predicts that Z goes to zero as the 
resonance is approached, since the size of the dressed 
molecules diverges at resonance and produces a vanishing 
overlap with the excited molecules. The measured quan- 
tity, however, continues smoothly through resonance, de- 
creasing exponentially with increasing field. Although 
the closed-channel fraction is finite and measurable, its 
magnitude above resonance is sufficiently small, <10 -5 , 
that the expectation of the number of closed-channel 
molecules is less than one. 
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FIG. 4: Z vs. B. The closed circles represent the value of 
Z extracted from measured values of Y. In the case of the 
920 G point, the loss is not exponential, but the value of V 
is taken to be the initial loss rate. The uncertainty in Z is 
approximately equal to the size of the closed circles, and is 
due mainly to uncertainty in the probe laser intensity. The 
dotted line shows a comparison with results obtained from a 
coupled channels calculation [l5j . The vertical dashed lines 
represent the boundaries of the strongly-interacting regime, 
fcjp|a| > 1, where &f is evaluated using typical values of N 
at the low and high field extremes. Although shot-to-shot 
variations in N are 30%, the average value of N at each field 
is between 13,000 and 90,000 due to day to day variations. 
Tp is between 200 and 600 nK due to differences in N as well 
as the trap frequencies. For all the data, T < T c and for the 
points above 850 G, T < 0.5T C , where T c refers to the critical 
temperature at 695 G. The gas is expected to be adiabatically 
cooled by the increasing field ramp for fields above 755 G 



The small closed channel fraction suggests compar- 
ison with a single channel model. We note that 
r is proportional to the local pair correlation func- 
tion G2(r, r) = (■0j(r)'!/'|(r)'0-f (r)-0|(r)), where $f and 

i{>l are the fermionic field operators for atoms in dif- 
ferent internal states. In the mean- field approxima- 
tion G2 may be factorized as C?2(r, r) = n 2 (r) + 
(t)) (r)tpi(r)), where the first term is the 
Hartree term with atom density n(r) = n-\{r) = ft^(r). 
This term gives rise to a slow two-body photoassociation 
process as was observed in the high temperature data of 
Fig. |21 The second term is non-zero only for correlated 
pairs and is proportional to |A| 2 , the square of the or- 



der parameter. In the BCS limit, |A| 2 cx e F e~ 7r/ ^ kF ^\ 
whereas in the BEC limit, |A| 2 oc ep/(kpa) 0, which 
is simply proportional to n(r), and produces a rapid one- 
body loss. In Fig. the data is compared with these 
functional forms. The fact that the data have the cor- 
rect dependence on (fcpa) -1 in the BEC and BCS limits 
is suggestive that such an approach has captured the es- 
sential physics. Data obtained by starting with higher 
initial temperatures do not fit the curves given in Fig. [S] 
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FIG. 5: Comparison of T with JA| 2 . T is expressed in units 
of fi 2 /7 making it equivalent to Z plotted in Fig. 01 The 
dashed lines correspond to evaluations of |A| 2 in the BCS 
and the BEC limits integrated over a Thomas-Fermi density 
profile. They have been scaled by the same factor to give the 
best fit on the BEC side. 

Note that non-condensed pairs will give rise to a sim- 
ilar factorization of Gi but it is not expected to have the 
same dependence on kpa. Although the order parameter 
presented in Fig.[S]comes from a single-channel model, it 
is the underlying closed-channel part, albeit small, that 
gives rise to the detectable signal, and indeed to the res- 
onance itself. 

We have shown that the singlet closed-channel com- 
ponent of the dressed molecules/pairs is quite small, less 
than 10~ 3 , throughout the strongly interacting regime. 
This result strongly supports the contention of universal- 
ity for broad resonances. Contrary to expectations from 
two-body physics, Z does not vanish for fields above reso- 
nance, but rather is continuous throughout the resonance 
regime. Furthermore, the molecular probe distinguishes 
between the paired and unpaired components of the gas 
and has enabled measurement of the order parameter in 
BEC and BCS limits, and throughout the strongly inter- 
acting regime. 
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